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Some results highlighting the status of a new version of a cluster fragmentation 
model for the Monte Carlo event generator Sherpa are presented. In its present 
version this model is capable of simulating e+e~ annihilation events into light- 
quark and gluon jets. We compare results for different multiplicity and momentum 
distributions to available SLD and LEP data as well as to results obtained with 
Herwig and Pythia. 



1 Introduction 

In view of the LHC starting to provide pp collisions at y/s — 14 TeV in 2007 
and a future e+e~ linear collider in the TeV range, new tools for the Monte 
Carlo simulation of multi-hadron final states are required. These tools will 
once be faced with the most precise measurements in high energy particle 
physics - this aspect dictating their improved physics content. In addition 
transparency and maintenance of these codes become an issue and therefore the 
object oriented language C-I--I- was chosen to write them. Beside the rewrites 
of the well-established tools Pythia ^ and Herwig ^, namely Pythia? ^ and 
Herwig H — 1-"^, another approach is available with the program Sherpa^. One 
of the striking features of Sherpa is the inclusion of the CKKW prescription 
to combine multi-jet matrix elements with parton showers ^ . This method 
allows a consistent description of multi-jet final states and a combination of 
such higher order calculations with the non-perturbative regime of hadron 
production in an universal manner. Even though the public Sherpa version ^ 
contains an interface to the Pythia string fragmentation, strong efforts have 
been made to develop a fragmentation model for Sherpa ^ that relies on the 
cluster fragmentation ansatz used within Herwig^. 



2 The Cluster Model 



A typical cluster fragmentation model consists of two parts: first, primary 
clusters are formed and, following this, such colour neutral states decay into 
hadrons and/or secondary clusters, which, in turn, decay further. 
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PYTHIA-e.Kttds) 


19.84 


16.72 


2.010 


0.856 


HERWIG-S.Ktids) 


18.86 


15.37 


1.693 


1.568 


SHERPACK 


20.15 


16.83 


2.018 


1.047 


DELPHI 


19.94 ± 0.34 


16.84 ± 0.87 


2.02 ± 0.07 


1.07 ± 0.05 


SLD 11 


20.048 ± 0.316 


16.579 ± 0.304 


2.000 ± 0.068 


1.094 ± 0.043 



Table 1: Overall mean charged-particlc multiplicity, and production rates of charged pions, 
charged kaons and (anti)protons in uds events at the Z'^-peak. 

In the first stop the partons emerging from tlie parton sliower are brouglit 
to tlieir constituent masses; tliis includes a finite gluon mass. The gluon tlien 
is forced to split into a liglit quark-antiquark (qq) or antidiquark-diquark (DD) 
pair. The resulting triplet and antitriplet states arc combined to obtain colour- 
neutral clusters. The model thereby allows for the incorporation of soft colour 
reconnection effects, which lead to configurations that are beyond the planar 
structure given by the parton shower evolution. However, these configurations 
are suppressed by a combined weight of and a kinematic function. After 

this first step four different cluster types can arise, mesonic {qq and DD), 
baryonic {qD)a,nd antibaryonic (Dq). 

In the second step the primary clusters, continuously distributed in mass 
with a peak at low cluster masses, have to be transformed into observable 
hadrons featuring a discrete mass spectrum. This is achieved by binary cluster 
decays and converting individual clusters into single primary hadrons. The 
model so far does not incorporate the subsequent decay of unstable hadrons. 
This task is still handled by the corresponding Pythia routines. Beside the 
request of locality and low momentum transfer in cluster decays, the model 
relies on a dynamic separation of clusters and hadrons. This implies that ac- 
cording to the flavour of its constituents a cluster is supposed to be a hadron, 
if its mass is below a certain threshold. Similar to the case of cluster forma- 
tion, the model for cluster decays incorporates the possibility for soft colour 
reconnection. Thereby a quark-antiquark pair or antidiquark-diquark pair pro- 
duced to disintegrate the cluster can recombine. Again these configurations 
are suppressed according to and a kinematic weight. 

The model so far is restricted to the fragmentation of light quarks (uds) 
and gluons produced in e+e" annihilation. Nevertheless various observables 
and distributions can be studied in order to validate the model. Figure 1 shows 
the mass distribution of the primary clusters at e+e~ collisions for different 
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energies proving the universality of the approach. Table 1 contains the mean 
multiplicities of tt^, and p, p in comparison with experimental uds results 
and Pythia and Herwig. In Figure 2 the charged particle hemisphere multi- 
plicity distribution and the charged particle scaled momentum prediction are 
compared to OPAL^^ and SLD^^ data, respectively, again taking into account 
uds events only. 

3 Conclusion 

The cluster-hadronization model developed for Shcrpa proved to work suc- 
cessfully for e+e^ annihilation events into light-quark and giuon jets. First 
tests show a satisfactory agreement with experimental data. Some cluster 
model shortcomings, such as the too low charged-particle multiplicities, could 
be cured; and the spectrum of the scaled momentum could be improved. The 
model will soon be extended to cover heavy-quark hadronization and the frag- 
mentation of beam remnants in hadron-hadron collisions. 
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Primary cluster mass distribution with CRM 
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Figure 1: The primary cluster mass distribution in e+e" annihilation events for different 
centre of mass energies using the colour reconnection model (CRM). 



Charged particle hemisphere multiplicity distribution Charged particle scaled momentum distribution 




Figure 2: Predictions for the hemisphere multiplicity distribution (left) and the scaled mo- 
mentum distribution (right) of charged particles considering the light-quark sector only. 
Results are compared to OPAL and SLD data, respectively. 
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